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技术 方法 
基于 密集 特征 匹配 的 胸 片 肺 野 目 动 分 割 


祭 广 南 , 陈 莹 筷 , 钟 丽 明 , 阳 维 , 冯 前 进 
南方 医科 大 学 生物 医学 工程 学 院 , 广 东 广州 510515 


摘要 :目的 X 线 胸 片 中 肺 野 的 准确 分 割 是 胸 片 图 像 自动 分 析 的 必要 步骤 。 本 文采 用 局 部 特征 的 密集 匹配 和 标号 融合 进行 胸 片 
肺 野 的 自动 分 割 。 方 法 对 于 输入 的 待 分 割 胸 片 ,基于 每 个 像素 点 提取 密集 SIFT 描 述 子 和 图 像 块 作为 局 部 特征 ,采用 密集 匹配 
直接 在 整个 参考 图 像 特征 集合 中 快速 搜索 近邻 ;密集 匹配 分 为 受 限 的 随机 初始 化 .近邻 场 传播 和 受 限 的 随机 搜索 三 步 ,并 数 次 
和 迭代 后 两 步 。 利 用 匹配 得 到 的 近邻 场 ,提取 标号 图 像 块 并 进行 加 权 , 权 重 为 匹配 的 相似 度 , 最 后 重组 为 肺 野 的 概率 图 ,经 闵 值 化 
处 理 即 可 得 到 肺 野 的 分 割 。 结 果 在 公开 的 JSRT 胸 片 图 像 数据 集 上 进行 测试 ,本 文 方法 的 Jaccard 指 标 可 达 95.5%。 结 论 利用 
局 部 特征 的 密集 匹配 和 标号 融合 能 取得 准确 性 高 且 稳 定 的 胸 片 肺 野 分 割 效 果 , 并 且 优 于 当前 最 好 的 胸 片 肺 野 分 割 方法 。 
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Automatic segmentation of lung fields in chest radiographs based on dense matching of 


local features 
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Abstract: Objective Accurate segmentation of lung fields in chest radiographs (CXR) is very useful for automatic analysis of 
CXR. In this work, we propose to use dense matching of local features and label fusion to automatically segment the lung 
fields in CXR. Methods For an input CXR, the dense Scale Invariant Feature Transform (SIFT) descriptors and raw image 
patches were extracted as the local features for each pixel. The nearest neighbors of the local features were then quickly 
searched by dense matching directly from the whole feature dataset of the reference images. The dense matching included 
three steps: limited random initialization, propagation of nearest neighbor field, and limited random search, with iteration of 
the last two steps for several times. The label image patches for each pixel were extracted according to the nearest neighbor 
field and weighted by the matching similarity. Finally, the weighted label patches were rearranged as the label class probability 
image of the input CXR, from which thresholds were obtained for segmentation of the lung fields. Results The Jaccard index of 
the proposed method reached 95.5% on the public JSRT dataset. Conclusion A high accuracy and robustness can be obtained 
by adopting dense matching of local features and label fusion to segment the lung fields in CXR, and the result is better than 
that of current segmentation method. 
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X 线 胸 片 是 目前 临床 上 广泛 使 用 的 医学 成 像 技 术 ， 
对 于 肺 部 疾病 的 分 析 与 诊断 具有 重要 的 临床 意义 。 对 
胸 片 中 肺 野 的 准确 分 割 是 对 疾病 的 分 析 与 诊断 的 重要 
前 提 。 目前, 随 着 计算 机 辅助 诊断 的 普及 ,已 经 有 很 多 
肺 野 自 动 分 割 算法 被 提出 来 了 ,分 割 准 确 度 与 速度 也 越 
来 越 高 ,但 仍 存在 一 些 难点 :锁骨 和 肋骨 窗 的 强 边界 , 容 
易 导 致 局 部 极 小 值 ; 极 小 的 胸 甩 和 角 难 以 准确 分 割 等 。 
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目前 胸 片 肺 时 分 制 算法 主要 有 基于 规则 的 分 割 趾 、 
像素 分 类 ,活动 形状 模型 等 。 基 于 规则 的 肺 时 分割 , 通 
篆 只 得 到 较为 粗糙 的 肺 野 轮 廓 ,可 用 做 活动 形状 模型 等 
的 初始 解 。 基 于 像素 分 类 的 分 割 取决 于 分 类 器 预测 性 
能 ,一 般 难 以 准确 定位 肺 野 边界 。 形 变 模 型 汪 由 于 具有 
灵活 的 形变 特性 而 广泛 应 用 于 医学 图 像 分 制 * ,但 也 存 
在 一 些 不 足 :(1) 对 初始 解 敏 感 , 易 陷入 局 部 极 小 值 点 ; 
(2) 对 算法 参数 敏感 ,造成 性 能 不 稳定 。 

为 实现 准确 可靠 的 肺 野 分 割 ,本 文 提 出 了 一 种 基 
于 图 像 特 征 密集 匹配 和 标号 融合 的 肺 野 分 割 方法 。 提 
取 待 分 割 胸 片 图 像 的 局 部 特征 ,将 像素 点 的 局 部 特征 与 
肺 野 已 知 的 参考 图 像 集中 的 局 部 特征 进行 匹配 ,确定 特 
征 在 参考 图 像 集 中 的 最 近邻 ,并 由 此 将 参考 图 像 中 的 肺 
野 标号 信息 通过 标号 融合 迁移 至 待 分 割 的 胸 片 图 像 ,以 
实现 肺 野 的 自动 分 制 |。 
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1 材料 和 方法 
1.1 材料 

实验 数据 采用 公开 的 JSRT 胸 片 图 像 数据 集 忆 ,和 包 
括 154 幅 肺 结 节 图 像 和 93 幅 无 结 节 图 像 , 其 中 每 幅 胸 片 
均 有 专家 手工 勾画 的 肺 野 轮 廓 。JSRT 数 据 集中 胸 片 原 
始 图 像 大 小 为 2048x2048 ,像素 大 小 0.175 mm, 灰 度 深 
度 为 12 比特 。 为 了 与 文献 结果 比较 和 减少 运算 时 间 ， 
原始 胸 片 图 像 降 采样 为 256x256 大 小 进行 实验 和 性 能 
评价 。 本 文中 算法 采用 Matlab 编 程 实现 。 
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图 1 本 文 方法 基本 流程 
Fig.1 Flow chart of the proposed method. 


1.2 方法 

1.2.1 方法 的 基本 流程 本 文 方法 分 为 两 部 分 实现 : 构 
建 参考 图 像 特征 集合 和 通过 密集 特征 匹配 分 割 肺 野 。 
算法 基本 流程 见 图 1。 给 定 参 考 图 像 集合 及 其 对 应 的 分 
割 标号 图 集合 ,参考 图 像 特征 集合 构建 过 程 如 下 :将 
所 有 参考 图 像 进 行 灰 度 归 一 化 处 理 ;@ 取 所 有 图 像 中 肺 
野外 接 和 矩形 内 的 区 域 ,将 其 长 宽 变 换 为 固定 之 后 合 加 平 
均 ,形成 模板 图 像 T;@) 将 所 有 参考 图 像 与 模板 图 像 T 进 
行 对 齐 , 对 参考 图 像 标 号 图 像 作 相应 的 空间 变换 ,提取 
局 部 特征 描述 子 并 降 维 , 得 到 参考 图 像 特 征集 合 以 及 对 
应 的 分 割 标号 图 像 集合 。 

在 构建 参考 图 像 特征 几何 完成 后 ,对 一 给 定 的 胸 片 
图 像 I, 对 其 进行 分 割 的 基本 过 程 如 下 :QD 对 输入 图 像 I 
进行 灰 度 归 一 化 处 理 ; 忆 将 图 像 1 与 模板 图 像 T 进 行 对 
齐 ;@@ 提 取 图 像 I1 的 局 部 描述 子 并 降 维 ;(@ 进行 K 次 密 
集 局 部 特征 匹配 ,得 到 开 个 近邻 场 ;@ 利 用 步骤 4 得 到 
的 近邻 场 和 局 部 特征 的 相似 性 进行 标号 融合 ;@) 将 标号 
融合 结果 进行 空间 逆 变 换 , 使 得 标号 融合 结果 与 原始 输 
和信 图像 在 空间 上 一 致 ;@ 将 标号 融合 结果 进行 闽 值 化 处 
理 得 到 肺 野 分 制 结果 , 阔 值 可 设 为 0.5。 

1.2.2 灰 度 归 一 化 和 空间 对 齐 为 减 小 胸 片 之 间 对 比 度 
差异 对 特征 匹配 造成 的 影响 ,对 胸 片 灰 度 进行 归 一 化 处 
理 。 设 胸 片 灰 度 直方 图 中 的 最 后 一 个 和 第 一 个 峰值 对 
应 的 灰 度 分 别 为 pi 和 ps, 归 一 化 处 理 后 的 图 像 灰 度 为 
Pi=(pi-Prin)/(Pmax -Pmin), 其 中 p 表 示 像 素 i 的 灰 度 值 。 


为 了 将 特征 匹配 的 搜索 范围 限制 在 较 小 区 域内 , 需 将 胸 
片 图 像 与 模板 进行 空间 对 齐 , 使 得 对 齐 后 的 胸 片 中 肺 野 
位 于 大 致 相同 的 位 置 。 本 文采 用 仿 射 变换 作为 形变 模 
型 ,通过 最 大 化 增强 相关 系数 “进行 胸 片 与 模板 的 对 齐 。 
1.2.3 图 像 块 匹配 图 像 块 匹配 是 指 对 图 像 块 在 图 像 中 
搜索 其 最 相似 的 图 像 块 (最 近邻 ) ,图 像 块 与 其 近邻 之 间 
的 位 移 关 系 由 最 近邻 场 表 达 。PatchMatch "是 一 种 高 效 
的 图 像 块 匹配 算法 。 由 于 图 像 的 一 致 性 , 相 邻 的 块 与 它 
们 的 最 近邻 可 能 具有 同样 的 位 移 向 量 ,因此 PatchMatch 
算法 可 以 将 最 近邻 场 传 播 到 相 邻 的 区 域 。PatchMatch 
算法 分 为 三 步 :随机 初始 化 .近邻 场 传播 .随机 搜索 。 随 
机 初始 化 指 的 是 对 图 像 A 中 的 每 一 个 图 像 块 都 在 图 像 
B 中 随机 地 指定 一 个 最 近邻 ,得 到 初始 最 近邻 场 ;近邻 
场 传播 可 以 把 好 的 位 移 关 系 传 到 周围 的 区 域 , 从 而 逐步 
优化 最 近邻 场 。 以 像素 点 psy 为 例 , 设 4i5 是 图 A 中 以 像 
素 点 paj 为 中 心 取 的 一 个 小 块 , 当 传 播 到 点 psy 时 , 设 块 
Bw 是 块 46y 的 最 近邻 , 且 与 465 相 邻 的 块 461) 4eio 都 
已 经 得 到 了 更 新 的 最 近邻 ,分 别 设 为 B Ba, 则 Bu 人 
8 有 可 能 是 4 的 最 近邻 ,分 别 计算 4 与 By Be 
Bu 这 三 个 块 的 欧式 距离 ,其 中 距离 最 小 的 则 作为 As 
传播 后 的 最 近邻 ,这 里 假设 Bi ,为 4 传播 后 的 最 近 
邻 ; 随 机 搜索 指 的 是 在 以 块 B6i ,的 中 心 像素 点 po.w 为 
中 心 的 搜索 窗 内 随机 选 出 一 个 块 B,, 然 后 再 从 中 择优 
选 出 一 个 作为 A j) 的 最 近邻 » 搜索 窗 的 大 小 依次 为 
W=mxnx(1/4) ,m 和 n 分 别 表示 图 像 4 的 高 . 宽 ,1 是 
随机 搜索 次 数 ,直到 W=1; 最 后 多 次 迭代 近邻 场 传播 和 
随机 搜索 的 这 两 步 ,最 近邻 场 会 逐渐 收敛 。 由 于 随机 初 
台 化 后 ,图 像 4 中 至 少 有 一 个 块 得 到 了 较 好 的 最 近邻 的 
这 种 可 能 性 很 大 ,所 以 ,经 过 多 次 迭代 的 传播 和 随机 搜 
索 , 最 后 4 中 绝 大 部 分 的 块 都 能 找到 合理 的 最 近邻 。 
PatchMatch 算 法 可 以 高 效 求 解 两 图 像 间 图 像 块 的 
近邻 场 。 但 当 参 考 图 像 数 量 较 多 时 ,逐一 将 输入 图 像 与 
参考 图 像 进行 块 匹配 是 不 实际 的 。 为 了 解决 这 一 问 
题 , Golland "提出 了 优化 的 块 匹 配 标号 融合 算法 ,用 于 
解放 结构 的 分 制 。 值 得 注意 的 是 , 块 号 配 算法 的 复杂 度 
取决 于 图 像 A 的 大 小 而 不 是 B, 优 化 的 块 匹配 算法 是 从 
参考 图 像 集中 搜索 个 近邻 ,而 不 是 只 从 与 图 像 A 相 似 
的 图 像 中 找 一 个 近邻 。 标 号 融合 的 时 候 , 每 个 近邻 的 
标号 的 权重 取决 于 匹配 对 之 间 的 相似 度 , 跟 块 匹 配 算 法 
一 样 ,不 是 只 取 中 心 像素 点 的 标号 ,而 是 基于 中 心 点 取 
了 一 个 块 。 另 外 ,随机 初始 化 和 随机 搜索 都 被 限制 于 一 
个 固定 大 小 的 搜索 窗 内 。 由 于 不 用 预先 检索 出 相似 图 
像 8, 这 个 优化 算法 效率 更 高 ,而 且 由 于 有 搜索 窗 的 限 
制 , 减 小 了 误 匹 配 率 , 实 现 了 更 高 的 准确 率 。 
1.2.4 特征 密集 匹配 与 标号 融合 对 于 胸 片 图 像 , 由 于 
图 像 采 集 条 件 和 病人 个 体 差 异 ,不 同 胸 片 肺 野 灰 度 和 对 
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比 度 差异 较 大 。 直 接 通过 图 像 块 匹配 和 标号 融合 ,难以 
得 到 可 靠 的 分 割 结果 。 为 了 解决 这 个 问题 ,本 文采 用 原 
始 图 像 块 和 密集 尺度 不 变 特 征 转 换 描述 子 (SIFT)" 作 
为 图 像 的 局 部 特征 。 原 始 像素 块 是 一 种 具有 较 强 区 分 
能 力 的 描述 子 , 但 对 噪声 和 形变 很 敏感 ,为 了 平衡 描述 
子 的 不 变性 和 区 分 能 力 , 原 始 图 像 块 的 尺寸 不 能 过 大 。 
SIFT 描 述 子 在 计算 机 视觉 有 着 广泛 的 应 用 ,具有 很 好 
的 旋转 不 变性 ,通常 用 于 描述 局 部 区 域 的 形状 特征 ,对 
于 图 像 的 形变 和 对 比 度 能 够 保证 一 定 的 不 变性 。 本 文 
中 ,原始 图 像 块 的 大 小 取 9x9,SIFT 描 述 子 的 维度 是 
128, 局 部 特征 的 维 数 为 209。 为 了 减 小 计算 复杂 度 , 采 
用 了 主 成 分 分 析 " 对 特征 向 量 进行 了 降 维 。 本 文中 ,局 
部 特征 的 维 数 降 为 40, 保 持 总 方差 的 95%。 

特征 匹配 与 块 匹配 不 同 ,对 于 输入 图 像 的 每 一 局 部 
特征 ,直接 在 所 有 参考 图 像 的 局 部 特征 中 搜索 其 对 应 的 
K 个 近邻 。 设 共有 NN 幅 参 考 图 像 , 参 考 图 像 大 小 为 mx 
n, 局 部 特征 的 维 数 为 DD, 则 参考 图 像 特征 集合 可 用 大 小 
为 DxNxmxn 的 四 阶 张 量 表达 ,参考 图 像 特征 集合 可 通 
过 离线 方式 提取 并 保存 。 特 征 匹配 搜索 近邻 的 步骤 与 
PatchMatch 算法 相似 ,分 3 步 进 行 迭 代 : 受 限 随机 初始 
化 .近邻 场 传播 . 受 限 随机 搜索 ,其 中 近邻 场 传播 仅 在 图 
像 平面 内 进行 。 

如 图 2 所 示 ,右边 的 大 长 方 体 表 示 特 征 向 量 集 , 圆 
点 代表 该 点 对 应 的 特征 , 浅 蓝 色 方 框 表示 固定 大 小 的 搜 
索 窗 ,表示 参考 图 像 序号 ,范围 是 1 到 N。 随 机 初始 化 
时 ,输入 图 像 特征 在 参考 图 像 特征 集合 中 随机 指定 一 个 
近邻 。 比 如 图 2 所 示 , 输 入 图 像 相 邻 点 的 近邻 可 以 位 于 
不 同 参考 图 像 中 。 近 邻 场 传 播 的 方法 与 PatchMatch 算 
法 类 似 , 仅 在 参考 图 像 的 图 像 平 面 内 进行 传播 , 即 近邻 
场 传 播 时 ,近邻 场 对 应 的 参考 图 像 序号 保持 不 变 。 随 机 
搜索 时 ,对 输入 图 像 任 一 点 随机 选择 参考 图 像 序号 ,在 
以 当前 近邻 为 中 心 的 搜索 窗 内 随机 选取 一 点 ,如果 随 机 
选取 的 特征 比 当 前 近邻 与 输入 图 像 特征 之 间 的 距离 小 ， 
则 将 随机 选取 的 近邻 作为 更 新 后 的 近邻 。 特 征 之 间距 
离 采 用 欧式 距离 。 

重复 特征 匹配 过 程 K 次 ,对 于 输入 图 像 可 以 得 到 KK 
个 近邻 场 (nn/) ,记录 了 输入 图 像 中 像素 点 特征 在 参考 
图 像 特征 集合 中 天 个 近邻 的 空间 位 置 和 参考 图 像 序号 
信息 。 在 标号 融合 过 程 中 ,本 文采 用 基于 图 像 块 的 标号 
融合 方法 "。 设 标号 图 像 中 像素 值 为 0 或 者 1, 其 中 1 
表示 像素 点 为 肺 野 。 对 于 输入 图 像 的 任 一 像素 点 (% 
y) ,在 标号 图 像 中 取 以 其 近邻 为 中 心 的 小 块 (如 5x5 大 
小 的 图 像 块 ), 并 以 点 (x, y) 的 特征 vy 与 其 近邻 的 特征 
ons 的 相似 度 作为 标号 块 的 权重 ,由 开 近 邻 加 权 得 到 
的 标号 块 为 


天 天 
lo,y) 二 > iD > mw: (1) 
i=1 i=1 


其 中 lonfi yy) 表 示 以 点 (x, 7) 在 第 i 个 近邻 场 nnf; 
对 应 的 近 人 Bnnfi oy 为 中 心 取 的 标号 块 ,权重 


0 oe 
w;=exp 


, SSDQV1,v2)= 
0 


yl -20) 家 示 两 个 D 维 向 量 的 欧式 距离 的 平方 , 另 
到 三 并 


外 ,s 的 取 值 要 适中 ,过 大 则 权重 都 约 为 1, 过 小 则 权重 都 
约 为 0, 实 验 中 可 取 为 待 分 割 图 像 中 所 有 像素 点 与 其 近 
邻 的 距离 的 均值 。 最 后 ,将 加 权 的 标号 块 进行 重 排 和 县 
加 为 输入 图 像 对 应 的 标号 图 ,并 进行 归 一 化 处 理 。 

为 了 得 到 与 原始 输入 胸 片 空间 位 置 对 应 的 分 割 ,还 
需 将 标号 融合 结果 用 模板 匹配 的 空间 变换 进行 逆 变 
换 。 标 号 融合 的 结果 可 以 看 做 肺 野 的 概率 图 ,对 其 进行 图 
值 化 或 使 用 图 割 方法 “进行 处 理 即 可 得 到 肺 野 的 分 割 。 


Features database 


Input image 


图 2 密集 特征 匹配 示意 
Fig.2 General view of dense matching of local features. Letft: 
Input Image, Right: Features Database of reference images. 


2 实验 和 结果 
2.1 实验 

实验 中 ,将 JSRT 胸 片 数 据 集 分 为 大 小 相同 的 两 个 
子 集 ,分 别 用 作 训 练 集 和 测试 集 。 肺 野 分 割 性 能 指标 包 
括 Jaccard 指 标 ( 开 ) 和 平均 边界 距离 (MBD),Jaccard 值 
越 大 .MBD 值 越 小 则 表明 分 割 的 性 能 越 好 ,其 中 : 


_IGTN SM| 
Jaccard = GTUSM| (2) 
MBD(GTB, SMB) = 
SD(gi,SMB) >)D(s,,GTB) 
1 i 了 
2 As 四 N, 


这 里 GT、SM 分 别 表 示 手 动 分 割 的 和 自动 分 割 的 
标号 图 ,1 为 肺 野 ,0 为 背景 ,公式 (2) 计 算 的 是 这 两 个 标 
号 图 的 肺 野 区 域 的 重 羡 率 , 重 羞 率 越 高 表示 分 割 | 效果 越 
好 ;公式 (3) 计 算 的 是 GT、SM 中 肺 时 的 边界 GTB、SMB 
之 间 的 欧式 距离 的 平均 值 , 当 手 动 分 割 的 肺 野 区 域 与 自 
动 分 割 的 完全 重 琶 ,Jaccard 值 为 1,MBD 值 为 0, 重 羡 率 
Jaccard 越 高 时 ,MBD 值 越 小 ,表示 分 割 效果 越 好 ， 
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图 3 近邻 数 不 同 的 标号 融合 结果 


Fig.3 Result of label fusion with different numbers of nearest neighbors. 


Dlg, SMB)= minllg; -s 咱 表示 GTB 上 的 任意 点 &; 到 


SMB 的 最 短 距离 ,N,N 分 别 表 示 GTB、SMB 上 像素 点 
的 总 个 数 。 
2.2 结果 

图 3 反映 了 特征 匹配 近邻 数 K 对 Jaccard 指 标的 影 


响 ,c-d 分 别 是 近邻 数 区 等 于 1.2.5.7.10 时 的 分 割 结 
果 。 网 4 显示 的 是 不 同 近 邻 数 时 标号 融合 后 的 结果 , 左 
中 至 右 近 邻 数 分 别 是 1.5.10。 可 以 看 出 近邻 数 增 加 时 ， 
分 割 得 到 的 肺 时 轮廓 更 平滑 ,分割 结果 更 准确 。 


入 


图 4 近邻 数 K 对 分 割 结果 的 影响 


Fig.4 Effect of K (number of nearest neighbors) on segmentation results. 


图 5$ 主要 反映 了 密集 匹配 过 程 中 迭代 次 数 对 分 割 
性 能 的 影响 , 横 轴 表示 近邻 数 。 图 中 可 见 Jaccard 值 在 
第 二 次 迭代 后 明显 提高 ,迭代 次 数 大 于 3 且 近 邻 数 大 于 
6 时 ,Jaccard 值 逐渐 趋 于 稳定 。 由 于 迭代 次 数 与 近邻 数 
越 大 ,需要 的 运算 时 间 会 增加 ,平衡 效率 和 分 割 性 能 ,本 
文选 取 的 迭代 次 数 为 3, 近 邻 数 为 10。 

图 6 反映 的 是 权重 系数 s 对 Jaccard 指 标的 影响 。 
当 s 很 小 ,比如 0.1 时 ,对 Jaccard 的 影响 较 大 ,尤其 是 近 
邻 数 较 小 的 时 候 ; 当 s 为 所 有 像素 点 与 其 近邻 的 欧式 距 
离 的 均值 ( 即 图 中 mean(Dist) , 约 为 9.5), 和 s 等 于 10 或 


100 时 ,Jaccard 值 几乎 完全 相等 ;所 以 ,为 省 去 计算 s 的 
时 间 , 本 文 直接 选取 s 为 10。 

图 7 是 JSRT 数 据 集中 部 分 图 像 自动 分 割 的 结果 与 
手工 分 割 结果 的 对 比 。 图 中 红线 和 绿 线 分 别 为 手工 勾 
画 和 本 文 方法 得 到 的 肺 野 轮廓 ;第 1 行 是 容易 分 割 日 效 
果 很 好 的 例子 ,第 2 行 是 较 难 分 割 且 分 制 效果 不 太 理 想 
的 例子 。 开 代表 Jaccard 指 标 。 表 1 列 出 了 本 文 方法 和 
其 他 肺 野 分 割 方法 在 JSRT 数 据 集 上 的 分 制 性 能 ,本 文 
方法 对 应 的 Jaccard 指 标 优 于 其 他 方法 。 表 2 列 出 了 各 
算法 用 相应 的 编程 语言 实现 时 运行 时 间 的 比较 。 
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23 456 78 9101112 图 6 权重 系数 的 影响 


Number'of nearestneighbors Fig.6 Effect of weight factor. 


图 5 密集 匹配 迭代 次 数 的 影响 


Fig.5 Effect of iterations of dense matching. 


于 有 
Ei” 


EF 
Ji=D.97043 


Fa 


JI=0.97433 Ji=097036 
有 
Jic0.95215 .66923 上 02 


图 7 JSRT 数 据 集中 的 图 像 分 割 结 果实 例 演 示 


Fig.7 Instances of segmentation results from JSRT database. 


为 了 验证 本 文 方法 的 有 效 性 和 推广 能 力 , 采 用 本 文 ”准确 。 昌 然 JSRT 数 据 集中 的 胸 片 都 是 通过 扫描 胶片 获 
方法 对 两 种 不 同 X 线 DR 设备 上 采集 的 数字 化 胸 片 进 ” 取 的 ,但 人 工分 割 的 结果 通过 本 文 方法 仍 可 很 好 地 迁移 
行 自动 分 割 。 图 8 显示 了 DR 胸 片 和 上 自动 分 割 肺 野 的 轮 ”至 DR 胸 片 的 分 拓 ,说 明 本 文 方法 具有 良好 的 推广 能 力 。 
廓 (红线 ), 可 以 看 出 分 割 结果 合理 ,并 且 肺 野 轮 廊 定 位 


3 讨论 


EY 


8 数字 胸 片 自动 分 割 实 例 


Fig.8 Instances of automatic segmentation in DR. 


结果 的 准确 性 。 从 实验 结果 看 ,本文 方法 在 一 定 程度 较 


胸 片 中 细小 的 胸 肪 角 和 受 锁骨 影响 的 肺 户 区 域 的 。 ”好 解决 了 这 些 问题 。 从 图 8 中 分 割 结果 较 差 的 图 像 来 
准确 分 割 比较 困难 ,并 且 肪 骨 徐 的 强 对 比 度 会 影响 分 割 。 看 ,严重 胀气 的 结肠 或 周围 组 织 被 错 分 为 肺 野 区 域 ,而 
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表 1 肺 野 分 割 算法 性 能 比较 


Tab.1 Comparison of performances of lung field segmentation 


algorithms 
Algorithms Jaccard (%) MBD (mm) 
Proposed method with Graph-cut 95.5+1.6 1 .3940.57 
post-processing 
Proposed method with threshold 95.441.6 1.4140.58 
post-processing 
Atlases with Non-rigid Registration"™ 95.4+1.5 1.3130.32 
Human observer'" 94.6+1.8 1.64+0.69 
下 and appearance Sparse 94.64+1.9 1.6740.76 
earning 
ASM tuned"" 92.7+3.2 2.30+1.03 


表 2 肺 野 分 割 运行 时 间 的 比较 


Tab.2 Comparison of execution time on lung fields segmentation 


Algorithms Execution time Programming language 
Proposed method with 12-148 Matlab 
Graph-cut post-processing 

Proposed method with 

threshold post-processing We Mebane 
ue Nonnad 20-258 C++ Matlab 
Registration 

Human observer'" - - 

Joint shape and appoatance 3528 Matlab 
sparse learning 

ASM tuned"®” ls C++ 


由 于 大 量 液体 存在 而 模糊 的 肺 野 区 域 被 错 分 为 非 肺 野 ， 
这 些 问 题 尚 需 进一步 解决 。 


4 结论 

本 文 提 出 了 一 种 采用 图 像 特征 密集 匹配 和 标号 融 
合 的 胸 片 肺 野 分 割 方法 。 与 一 般 的 两 图 像 间 块 匹配 方 
法 不 同 , 本 文采 用 密集 SIFT 和 原始 图 像 块 作为 图 像 的 
局 部 特征 ,直接 将 输入 图 像 特 征 与 参考 图 像 集合 的 特征 
进行 匹配 ,避免 了 参考 图 像 子 集 的 选择 问题 ,有 效 利 用 
所 有 参考 图 像 及 其 标号 的 信息 。 在 JSRT 胸 片 图 像 数据 
集 上 ,本 文 方法 得 到 了 目前 为 止 最 高 的 分 割 精 度 。 此 
外 ,实验 结果 还 表明 本 文 方法 具有 很 好 的 推广 能 力 。 在 
进一步 的 工作 中 ,我 们 将 对 特征 匹配 方法 进行 推广 应 
用 ,如 用 于 胸 片 中 结 节 的 自动 检测 和 肋骨 的 抑制 。 
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